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Abstract 

The statistical properties of random analytic functions ip(z) are investigated 
as a phase-space model for eigenfunctions of fully chaotic systems. We gen- 
eralize to the plane and to the hyperbolic plane a theorem concerning the 
equidistribution of the zeros of ip{z) previously demonstrated for a spherical 
phase space (SU(2) polynomials)). For systems with time reversal symme- 
try, the number of real roots is computed for the three geometries. In the 
semiclassical regime, the local correlation functions are shown to be univer- 
sal, independent of the system considered or the geometry of phase space. 
In particular, the autocorrelation function of ip is given by a Gaussian func- 
tion. The connections between this model and the Gaussian random function 

hypothesis as well as the random matrix theory are discussed. 
PACS numbers: 05.45. +b; 03.65.Ge; 03.65.Sq 
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I. INTRODUCTION 



There are two basic descriptions of the eigenf unctions of fully chaotic systems. On 
the one hand, there is the random matrix theory of Wigner and Dyson as a model of the 
quantum behavior of chaotic motion 0,0. In this model the Hilbert space is isotropic, i.e. 
an arbitrary eigenstate <p of the random matrix can point out in any direction of Hilbert 
space with equal probability. The only constraint comes from the normalization. Denoting 
by a = (ao, ax, ■ ■ ■ , a^) the components of (p in a given (and arbitrary) basis, the joint 
probability distribution for the amplitudes is M 
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where \S n \ = 2n n ^ 2 /Y{n/2) is the surface of a (n — l)-dimensional sphere of unit radius. We 
assume for the moment that a is a complex vector, this being the usual situation when the 
physical system has no time-reversal symmetry. By construction this theory concentrates 
on statistical aspects of the eigenstates, and the "ergodicity" in Hilbert space expressed by 
Eq. ([!]) is a direct consequence of the invariance of the ensemble of matrices under unitary 
transformations. 

On the other hand, a semiclassical theory has also been developed 0-0]. Unlike the 
random matrix theory, its aim is to provide a system-dependent description of the eigenstates 
based on the solutions of the classical dynamics. Adopting a phase-space representation of 
quantum mechanics, in the semiclassical limit % — > the Wigner transform W^(x, p, E) of 
an eigenstate ip of an ergodic system at energy E takes to leading order the simple form 

W4x,p,E)^Ar5[E-H(x,p)} , (2) 

where M is a normalization constant, (x,p) parametrizes the phase space and H(x,p) is the 
classical Hamiltonian. Because in the semiclassical limit the Wigner function presents wild 
oscillations, Eq.(^j) has to be interpreted as a smoothed version over a phase space region of 
the order of one Planck volume. 
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Besides the energy-shell ergodicity, Eq.(0) reveals a strong localization effect in phase 
space of which there is no trace in Eq.(|l]). Because of energy conservation is, in the 
semiclassical limit, exponentially concentrated on the energy shell. This means that when 
expanding an eigenstate of H in an arbitrary basis m of Hilbert space it will have com- 
ponents different from zero only on those states whose Wigner transform intersects - or 
is included into - the energy shell 8 [E — H(x,p)}, while the amplitudes will vanish expo- 
nentially fast otherwise. Therefore, in general for a given chaotic Hamiltonian system the 
expected distribution of the eigenstates's amplitudes is quite different from the isotropic 
rule (P, because random matrix theory ignores the phase-space restrictions imposed by en- 
ergy conservation. Eq.(|]) has therefore to be interpreted as a valid statement only for the 
amplitudes a m whose associated vector m intersects the energy shell. 

Having this in mind, our purpose is to further develop the connections between the two 
above mentioned descriptions. In particular, assuming a two-dimensional phase space we 
study the ergodicity of the eigenstates as implied by Eq.(|l|), as well as some correlation 
functions computed inside the energy shell. For that purpose, we use the Bargmann (or 
coherent state) representation which associates to each quantum state an analytic function 
if^(z) in the complex plane. This complex plane is interpreted as a parametrization of 
the energy shell manifold, and the three simplest one-dimensional complex manifolds are 
considered: the sphere (curvature +1), the plane (zero-curvature) and the hyperbolic plane 
or pseudosphere (curvature —1). After introducing the model in section II, we show in 
section III that the associated Bargmann functions given by Eq.(Q) below have, when the 
coefficients a m are distributed according to Eq.(|l|), a uniform distribution of zeros over the 
sphere, the plane and the hyperbolic plane, respectively. The uniformity of the distribution 
was shown in ||[| for the spherical geometry, and the results of this paper constitute a 
generalization to the plane and to the hyperbolic plane. Some related theorems extending 
some of the results to higher dimensional phase spaces were recently obtained independently 
in Ref. HTJ. 



We also consider the case when the W\ and SU(1, 1) series in Eq.(H) are truncated to a 



finite number of terms (thus introducing W\ and SU(1, 1) polynomials). It is shown that 
for those polynomials the distribution of roots is now uniform over a finite region (a disk) 
of the complex plane and the hyperbolic plane, respectively 

The equidistribution of roots holds only if the vector a is complex. When it is real (which 
corresponds to time- reversal symmetric systems), a concentration of roots occurs over the 
real axis in the complex plane. We show that for the three random functions considered 
the number of real roots located inside a disk centered at the origin containing N zeros is 
proportional to y/~N. 

Section IV is devoted to the study of the correlation properties of the random analytic 
functions between different points of the energy shell. We first argue that in the semiclassical 
limit all the correlation functions are locally universal, independent of the geometry of phase 
space. We then compute as a particular example the autocorrelation of ip(z). 

II. RANDOM ANALYTIC FUNCTIONS 

In order to study the statistical properties of eigenstates of chaotic systems, we first 
define an energy-shell (phase space) associated to the quantum mechanical system. We will 
not consider a specific dynamical system but make general statements which are expected 
to be valid for any system having at the classical level a strongly chaotic behavior. 

Concerning the energy shell, the minimum dimensionality required in order to have 
chaotic behavior at the classical level is two. This two-dimensional manifold M could be for 
example the Poincare surface of section of the energy shell of a two-freedoms Hamiltonian 
system. But we can alternatively choose a two-dimensional manifold and directly define 
on it a fully chaotic area-preserving map, thereby avoiding the problem of the reduction 
to the Poincare section. In both cases one can define a unitary operator which quantizes 
the classical flow. In the following we study properties of the eigenstates of this unitary 
operator. Because in principle it corresponds to a section of it, we will still call the manifold 
M. the "energy shell". 
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Independently of the definition of the dynamics, the three simplest possibilities for Ai are, 
in two dimensions, the sphere, the plane and the hyperbolic plane or pseudosphere JTT] (the 
cylinder and the torus being just periodicized versions of the plane), having curvature +1, 
and —1, respectively. The choice of the manifold will depend on the dynamical symmetry 
group underlying the system Hamiltonian under study. The three spaces considered are 
associated with the SU(2), W\ and SU(1, 1) symmetry, respectively. For example, for the 
motion of a spin in a variable magnetic field the appropriate phase-space manifold is the 



sphere (see Ref. |12|] for illustrations in different geometries). Area preserving chaotic maps 
have been studied in the plane and its periodicized versions (e.g. the cat, standard and 
kicked Harper maps) and the sphere (the kicked top) but not, to our knowledge, in the 
pseudosphere. 

The next step is to define a representation of the quantum states over the manifold Ai. 
This can be achieved for the three geometries under consideration by means of a phase-space 
coherent state (Bargmann) representation of quantum mechanics consisting of overcomplete 
and non-orthogonal systems of Hilbert space vectors. An arbitrary quantum state can be 



written in the coherent state basis as [ P~3"| , |T2" 
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where C l k = k\/l\(k — l)\ are Newton's binomial coefficients. For SU(2) we have iV = 2j, 
where j = 1/2, 1, 3/2, ■ ■ • is the usual spin label of the (2j + l)-dimensional representation. 
For SU(1, 1), L = 2k, with k — 1, 3/2, 2, • ■ ■ is the usual label of the discrete representation 
of the group [|T2|]. Unlike the case of W\ and SU(1, 1), the representation of the quantum 
state for SU(2) is finite dimensional, reflecting the fact that the group is compact. 

Eq. (HI) provides a realization of the Hilbert space in terms of analytic functions ip( z )- 
The complex variable appearing in Eq.(|3D spans in each case the associated manifold Ai. It 
is defined as 



cot (0/2) exp (iif) , 

(X - ip)/y/2 , 

tanh(r/2) exp (iip) 



SU(2) 
W l 

517(1,1) . 



(4) 



Here (#, </?), (x,p) and (r, </?) are the usual variables labeling these manifolds (see Fig.l). The 
measure of these three manifolds projected onto the z-plane is 

r i d 2 z 
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. 7T (1 — Zz) 2 

The notation d 2 uo where u is a complex variable stands for d(Reu)d(lmu) . 

Now comes the statistical ingredient entering in the definition of the model. If ip is an 
eigenstate of the quantization of the map acting on M. then, for a fully chaotic dynamics, 
we assume that the amplitudes {a m } in Eq.(|) are distributed according to Eq. ([[]). Of 
course this is an hypothesis whose validity has to be compared with the behavior of real 
dynamical systems. Including this element, we obtain therefore a phase-space model for 
chaotic eigenstates in terms of random analytic functions. For SU(2) polynomials this 



model was introduced in Ref. and further considered in [|T^,|T^,[|,|l^,[T7]] . 

Being analytic, the average properties of the random functions ij){z) cannot be stationary 
over the complex plane (unless the function is trivially a constant). For example, for 577(2) 
the average square modulus of ip(z) grows like 

N 
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(6) 
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The brackets indicate average over a uniform distribution for the coefficients, (f(a)) = 
J d 2 ao . . . d 2 aNV(a)f(a). For convenience we use for the computation of the averages a 
different normalization from that of Eq. (|l|) 
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We therefore have 



(if>{z)if>{z)) 



(flnfidji) &mn j (8) 

from which we obtained Eq.(^). This change of normalization has no incidence in the results 
reported in this paper since the quantities considered (cf Eqs.([l|) and (pOD ) are independent 
of the global normalization of i/j. 

An analogous computation can be done for the other two geometries, the plane and 
the pseudosphere. When computing the averages for these two geometries, we truncate the 
series in Eq.([|) to a finite value m = N and then we let N — > oo. We get 

(l + zz) N SU(2) 

exp (zz/h) Wi (9) 
(l-zz)- L SU(1,1) . 
As was pointed out in Ref. []18 |, the statistical properties of the zeros of ip{z) are un- 
affected if instead of the distribution (|7|) a Gaussian independent distribution for the coef- 
ficients is used. This is also true in the large N limit for any statistical property of i/j(z), 
because the distribution (|7]) also implies in that limit a Gaussian law for the distribution 
of a given coefficient. For most purposes, Eq.(|7]) (or alternatively ([!])) is equivalent to a 
Gaussian independent distribution. 

The random analytic model then implies that at each point of phase space the distribution 
P{ip) of the function ip{z) is Gaussian, with a local variance given precisely by Eq.(^) (in 
agreement with the old conjecture of Berry 0) 

Pty)= 1 1 2 exp(-|^)| 2 /2(|^)| 2 )) . (10) 
The relationship between the statistical properties of ip in different representations was 



recently considered in 19 . 



III. ERGODICITY 



The simplest question concerning the statistical properties of the functions ip(z) is related 
to the ergodicity. We would like to establish more clearly the connection between the Hilbert 
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space ergodicity of Eq.(P on the one hand and the phase-space ergodicity over the energy 
shell on the other. From the results of the previous section this connection is now easy to 
establish. 

As we mention, the functions ip{z) are not stationary in the complex z-plane. However, 
the quantity which is physically associated to the presence of the particle in a neighborhood 
of a phase space cell (and therefore more closely related to the notion of ergodicity) is not 
ij)(z) but rather the (quasi) probability density 

usually called the Husimi function. It is simply obtained from the Wigner transform through 
a Gaussian convolution. Exactly as for the average Wigner function Eq.(|2]), Eq. (|TT|) implies 
that the average of the Husimi function is constant in the complex z-plane. Because of this, 
the (quasi) probability TL^d^ is going to be stationary on the manifold M. (and not in the 
complex z-plane where the measure is not uniform (cf Eq.(|5])). 

In particular, the previous statement applies to the distribution of the zeros of H^, i.e. 
the density of zeros of TC^ should be uniform over the manifold Ai . This is going to be true 
for the zeros of ip{z) as well, because both functions have the same set of zeros [2D|. In fact, 



the uniformity over M. of the density of zeros p{£) = 5[ip(z)} \dip(z)/dz\ 2 of ip{z) can be 
directly checked for the three geometries from the general formula 



1 d 



2 



<pW> = -^MV>(*M*)>] • (12) 



From @ and ([12]) we find, for each geometry 

( N d 2 z N 



sin9 d9d<p SU(2) 



7T (1 + ZZ) 2 AlT 

(p(z)) d 2 z = | d 2 z/(nh) Wi (13) 

~i — ~ Z -\o = — sinh r drdtp SU(1, 1) . 

k 7T (1 — ZZ) 4:71 

For the sphere and the pseudosphere, the density has been transformed to the usual 
coordinates spanning the corresponding manifold. Eq.flTBD shows that the density of zeros 
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of the random--?/' (z) model is indeed uniform over the associated manifold Ai (this being 
true for arbitrary values of the parameters N, L or h). 

The uniformity of the density of zeros was also demonstrated for random analytic func- 
tions defined on the two-dimensional torus There exist also some related results ob- 



tained without using the random hypothesis. For example it was proved [|2T| that the 
Schnirelman property for the Husimi density implies a (weak) convergence of the density of 
zeros towards the uniform measure. The uniformity of the distribution of zeros was checked 
numerically for certain maps on the sphere |||9| and also for two-dimensional chaotic systems 



with W\ symmetry |22| 



In spite of its aesthetical beauty, Eq.(pT3[) has an unpleasant feature for the W\ and 
the »S77(1, 1) groups, since a uniform density of zeros for ip(z) over the whole (non-compact) 
manifold M. is not realistic from a dynamical point of view. For chaotic Hamiltonian systems, 
non-compact energy shells like the plane or the pseudosphere are somewhat artificial because 
energy conservation implies in general a compact manifold (or Poincare surface of section). 
A more appropriate model for the eigenfunctions of a system having a W\ or SU (1,1) 
symmetry would then be a truncated version of Eq.^) to a finite degree N. For the plane 
for example, we have the W\ polynomial 

m=0 Vmlh 

In practice the integer N at which the truncation is made fixes the approximate number 
of zeros lying in the classically allowed region of phase space at a fixed energy (outside 
which the associated Husimi function has a rapid decay). It is given by the semiclassical 
rule N m (area of classically allowed region) /27rh. The average density of roots for this 
polynomial can be computed from (12), with the result 



N 



where g(x) = (x /TV!)/ J2 (x m /m\). The function ( JIB] ) has a step-like shape with a constant 



m=0 



density of zeros 1 / (irh) inside a circle of radius ~ V Nh and a vanishing density outside (see 



Fig. 2). It tends, in the large A-limit, to the unit step function 0(1 — r) where r = \z\j\J Nh. 
The convergence is however slow since for large values of r we find a power law dependence 
(p(r)) « l/(7chNr A ) + C(l/r 6 ). 

We therefore recover for the W\ random polynomials the appealing feature that the 
Bargmann transform of a chaotic eigenstate has a uniform density of zeros in a finite region 
of the phase space plane (over which the associated Husimi function is stationary). The 
decay of the density outside the classically allowed region is also of interest, because it 
describes the evanescence of the eigenfunction ip in the classically forbidden region |23| . 



But we don't expect the decay predicted by Eq. ([l5|) to be very accurate due to the abrupt 
truncation of ip at m = N. 

Analogous results are obtained for the SU(1, 1) polynomials defined, like the W\ poly- 
nomial, by the truncation of the series in Eq.(^) to a finite value of m 

N 



m=0 

In the limit iV > L we find a constant density of zeros in the hyperbolic plane inside a 
region of size coshr ~ N/L. 

Time reversal symmetry 

The previous results are valid when the coefficients of the random functions are complex. 
This generally corresponds to systems without time reversal symmetry. For systems having 
that symmetry, the coefficients {a m } can be chosen to be real, and ip(z) satisfies now the 
functional equation 



ip(z) = V(z) • 

As discussed in detail in ||, the existence of a functional equation generically produces 
a concentration of a certain number Nr of roots over the symmetry line (here the real 
axis in the complex z-plane, labeled by x). The distribution and correlations of roots of 
real polynomials was studied in [J24| . For SU(2) polynomials, explicit expressions for the 
correlations existing between real roots where obtained in [JE 
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The number Nr = J(p(x)) dx of real roots of a polynomial whose coefficients are real 
random Gaussian independently distributed with variance <r^ may be computed from a 
general formula obtained by Kac. The average density of real roots is given by (18| 



1 JA(x)C(x) - B 2 (x) 



with 

N 



_2 2m 



A(s) = £ 

m=0 

B(x) = A'(x)/2 ( 18 ) 
C(x) = [A"(x)+A'(x)/x}/4: , 
where the prime indicates derivative with respect to x. 

For SU(2) polynomials we have cx^ = C™, = (1 + a; 2 )^ and from (TL7|) we obtain 



(p{x)) dx = By a stereographic projection x = cot (0/2) the density simplifies to 

(p(6)) d9 = ^-d9. After integration we obtain N R SU ^ = \Z~N, a result which was obtained 
in 0. 

For real W\ random functions, = l/(k\h ), A(x) = exp(x 2 /h) and (p(x)} dx = 
dx/(ir\/2h) (remember that the measure is Since for the complex W\ function we 



have, according to Eq . fll3|) , N roots in a disk of radius r = y2hN, we find that the number 

of real roots when the coefficients are real is N^ l] = ^^N. 

Finally, for SU(1, 1) we have a 2 m = C^ m _ x , A{x) = 1/(1 -x 2 ) N and (p(x)} dx = f 

Going to the hyperbolic plane via x = tanh(r/2), the density is simply written (p(r)) dr = 

^-dr. Since for the complex S77(l,l) function we have, according to Eq. (|T3|) , a number 

L of roots in a region of size r = arccosh(3), we find a number of real roots given by 
jy(SU(i,i)) _ arccosh(3) ^/^y 



IV. PHASE-SPACE CORRELATIONS UNIVERSALITY 



Besides the average density of roots, other relevant statistical properties are correla- 
tion functions of ip{z) or of its zeros. One might think that this quantities are geometry- 
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dependent. However, in the large- N limit, the three different geometries treated here tend 
to the same flat-space result. Indeed, for fixed m we have 

X m 

lim = hm C^ m _ x = — , (19) 

X^oo X^oo Tfl\ 

where X = N for the SU(2) group and X = L for the SU(1, 1). This implies that, with 
the scaling z — > yXz, the SU{2) and truncated £77(1, 1) polynomials tend to the truncated 
W\ polynomial in a region around the origin of the z-plane (the argument can be repeated 
for an arbitrary point lying inside the const ant- density region). This is easy to understand 
since in the large-iV limit the local density of zeros of the rescaled polynomials tend to 
diverge (the mean spacing between zeros in rescaled units goes like 1/y/N). On a scale 
of distances of a few mean spacings between zeros the space will look, locally, flat. As a 
consequence the local fluctuation properties of all these functions are, asymptotically, the 
same J25j. This can be rephrased as a general statement concerning the local correlation 
properties of eigenstates of fully chaotic systems (for points lying inside the energy shell): 
they are universal, independent of the system under consideration or the geometry of phase 
space. 

The n-point correlation functions among zeros have been computed analytically by J. 
Hannay for random 577(2) polynomials |L5]. In particular, it was found that there is a cubic 



repulsion between zeros. The previous considerations imply that in the large- N limit these 
correlation functions should be valid locally for any geometry. The two-point correlation 
function as well as a computation of the nearest neighbor spacing distribution for the random 
analytic model have been compared with numerical results obtained for several chaotic maps 



in different phase space geometries in [plj. Very good agreement has been found, thus 
supporting the universality of the correlations. Further support comes from a computation 
of the two-point correlation function and the nearest neighbor spacing distribution for the 



zeros of a two-dimensional chaotic system having a W\ symmetry P2" 

Concerning ip{z) itself (and not its zeros), the simplest quantity of interest is the auto- 
correlation function 
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As before, this quantity will be asymptotically independent of the geometry. For simplicity, 
we restrict the computation to the W\ case. Using Eqs.@, (||) and (||) one finds 

C(Z,z) = exp[-\Z\ 2 /2 + iIm(zO] . (21) 

Unlike the phase, the modulus of the autocorrelation is independent of z and is isotropic. 
This result should be compared with an analogous computation in configuration space [Q, 
where it was found that C is given by a Bessel function (we find that Eq. ( |2"ID holds also for 
systems with time reversal symmetry; analogous quantities for systems with and without 
time-reversal symmetry have been considered more recently in [p6fl ). 

V. CONCLUSIONS 

The distribution of the zeros of W\ random polynomials of degree N is strongly reminis- 
cent to that of the eigenvalues of a Gaussian ensemble of N x N complex matrices discussed 
by Ginibre (see chapter 15 in |T|]). Indeed, for large values of N both distributions have a 
uniform density in the complex plane inside a disk of radius \z\pd y/~N, and rapidly decaying 
to zero outside of it (for simplicity, we set h — 1 in Eq. (|T6|) ) . The distributions are however 
different. The density has a Gaussian tail for \z\ ^> in the case of the eigenvalues of 
random complex matrices, while it has a power law decay for W\ polynomials. Moreover, in 
the case of Ginibre the distribution may be interpreted as a two-dimensional Coulomb gas, 
while the two-dimensional gas associated to zeros of random polynomials includes, apart 



from the Coulomb interaction, iV-body terms |16[ . A detailed comparison between these 
two interacting systems was recently made in |27 . 



Several extensions of the results reported here are conceivable, in particular to higher 
dimensional complex manifolds. The random matrix distribution (P is in fact valid indepen- 
dently of the number of degrees of freedom of the underlying physical system. It is therefore 
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natural to expect that results analogous to those obtained here for the distribution of the 
zeros in two dimensional phase spaces also hold in higher dimensions (like for example for 
polynomials with several complex variables). And indeed a generalization to SU(m + 1) 
polynomials with m > 1 was recently done in Ref. [ |T0|j . 

Acknowledgments: I am indebted to O. Bohigas and A. Voros for many stimulating 
discussions. 
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Captions to Figures 



Figure 1: (a) The sphere and, (b) the pseudosphere (upper sheet of the hyperboloid 
— Xq + x\ + x\ = —1) having curvature +1 and —1, respectively. For SU(2), the definition 
of z corresponds to a stereographic projection which maps the unit sphere labelled by (9, ip) 
onto the equatorial complex plane labelled by z, the north pole of the sphere being the center 
of projection. In the case of SU(1, 1), a projection of the pseudosphere onto the unit disk 
in the complex z-plane located at x = is done through the point (x\, x 2 , x ) = (0, 0, — 1). 



Figure 2: The normalized density of zeros Eq.(p~5D for the W\ random polynomials as a 



function of r 




for N = 100. 
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